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ABSTRACT 

Correlations between the star formation rates (SFRs) of nearby galaxies (so- 
called galactic conformity) have been observed for projected separations up to 
4 Mpc, an effect not predicted by current semi-analytic models. We investigate 
correlations between the mass accretion rates (dM v i r /dt) of nearby halos as a 
potential physical origin for this effect. We find that pairs of host halos “know 
about” each others’ assembly histories even when their present-day separation is 
greater than thirty times the virial radius of either halo. These distances are far 
too large for direct interaction between the halos to explain the correlation in 
their dM vir /dt. Instead, halo pairs at these distances reside in the same large- 
scale tidal environment, which regulates dM v i r /dt for both halos. Larger halos are 
less affected by external forces, which naturally gives rise to a mass dependence of 
the halo conformity signal. SDSS measurements of galactic conformity exhibit a 
qualitatively similar dependence on stellar mass, including how the signal varies 
with distance. Based on the expectation that halo accretion and galaxy SFR are 
correlated, we predict the scale-, mass- and redshift-dependence of large-scale 
galactic conformity, finding that the signal should drop to undetectable levels by 
z > 1. These predictions are testable with current surveys to z ~ 1; confirmation 
would establish a strong correlation between dark matter halo accretion rate and 
central galaxy SFR. 

Key words: cosmology: theory — dark matter — galaxies: halos — galaxies: 
evolution — large-scale structure of universe 


1 INTRODUCTION 

In the ACDM paradigm, galaxies form at the centers of 
collapsed, gravitationally self-bound halos of dark mat¬ 
ter. Across most of cosmic time, galaxy stellar mass ex¬ 
hibits a tight statistical connection with dark matter halo 
mass (e.g., Conroy & Wechsler 2009; Moster et al. 2010, 
2013; Behroozi et al. 2010, 2013a,b; Firmani & Avila- 
Reese 2010; Leauthaud et al. 2012; Leitner 2012; Yang 
et al. 2012; Bethermin et al. 2013; Wang et al. 2013b; 
Mutch et al. 2013; Lu et al. 2014b). Yet, individual galax¬ 
ies enjoy a diversity of star formation rates (e.g., Brinch- 
mann et al. 2004; Noeske et al. 2007; Salim et al. 2007; 
Brammer et al. 2011; Moustakas et al. 2013; Muzzin et al. 
2013), which has proven more difficult to physically and 
self-consistently link to halo properties (Lu et al. 2014a; 
Genel et al. 2014; Behroozi et al. 2015). 

A recent insight into this link is galactic conformity : 
the tendency of neighboring galaxies to have similar spe¬ 
cific star formation rates (SSFRs), colors, gas fractions, 
and morphologies. This effect was first observed for satel¬ 


lites of larger (i.e., “central”) galaxies (Weinmann et al. 
2006; Wang et al. 2010; Kauffmann et al. 2010; Wang & 
White 2012; Robotham et al. 2013; Phillips et al. 2014; 
Knobel et al. 2015), but later also found to occur for 
galaxies separated by up to 4 Mpc in projection (Kauff¬ 
mann et al. 2013). Naively, this larger scale “2-halo con¬ 
formity” (see Hearin et al. 2014b) would be difficult to 
explain via internal baryonic processes alone, as these 
correlations extend over 1000 times the half-light radii of 
the galaxies and over 10 times the virial radii of the halos 
concerned. 1 Two-halo conformity is evidently quite chal¬ 
lenging to correctly predict: the signal in contemporary 
semi-analytic models (Guo et al. 2011) and the state-of- 
the-art Illustris hydrodynamical simulation (Bray et al., 


1 Observations of “1-halo conformity” also remain to be suc¬ 
cessfully explained, though this paper will chiefly be concerned 
with correlations beyond the halo’s virial radius. See, however, 
§4.3. 
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private communication) is nowhere near as strong as in 
SDSS observations. 

Although two-halo conformity in galaxy SSFR may 
seem puzzling, in fact many properties of dark matter 
halos are correlated on these scales (Hahn et al. 2007a,b; 
Wang et al. 2011; Shi et al. 2015). Such correlations are 
not at all puzzling, as nearby halos form out of the same 
over-dense patch in the cosmic density field. This phe¬ 
nomenon has been investigated extensively in the context 
of “assembly bias”—a dependence of the clustering of ha¬ 
los on properties besides mass (Gao et al. 2005; Wechsler 
et al. 2006; Wetzel et al. 2007; Wang et al. 2007; Kesel- 
man k Nusser 2007; Gao k White 2007; Croton et al. 
2007; Wu et al. 2008; Li et al. 2008; Dalai et al. 2008; 
Faltenbacher k White 2010; Lacerna k Padilla 2011, 
2012; Wang et al. 2013a; Zentner et al. 2014). Previous 
studies have also established long-range correlations be¬ 
tween tidal forces and internal halo properties and/or 
halo assembly history (Wang et al. 2007; Hahn et al. 
2007a,b, 2009; Wang et al. 2011; Behroozi et al. 2014; Shi 
et al. 2015). There is no shortage of halo properties with 
large-scale correlations. From this perspective, two-halo 
conformity—especially conformity between the SSFRs of 
central galaxies—should naturally arise if star-format ion 
history were coupled sufficiently strongly to one or more 
halo properties that exhibit large-scale correlations. 

The dark matter accretion rate of a halo is tightly 
coupled with the accretion rate of baryonic gas (see Wet¬ 
zel k Nagai 2014, for a recent demonstration). This tight 
coupling provides a theoretical basis for connecting the 
SFR of a galaxy to its halo’s dark matter accretion rate 
(dM V i r /dt). In this paper, we investigate whether a con¬ 
nection of this type could give rise to an SDSS-like con¬ 
formity signal, and predict how parameters of the confor¬ 
mity measurement (distance, mass, and redshift) affect 
the strength of the signal. 

This paper is organized as follows. We describe the 
N-body simulation we use in §2, and in §3 we demon¬ 
strate that halo accretion rates are correlated over the 
same range of scales as SDSS galaxy star-formation rates. 
We identify the root cause of this phenomenon as corre¬ 
lations in the tidal forces experienced by nearby halos. 
This results in a unique prediction for the halo mass and 
redshift-dependence of two-halo conformity under our hy¬ 
pothesis that galaxy star formation rate is coupled to halo 
accretion rate. We discuss tests of these predictions and 
implications in §4, and summarize our primary conclu¬ 
sions in §5. 


2 SIMULATION k HALOS 

We use the Bolshoi simulation (Klypin et al. 2011) for 
all our analysis. Bolshoi follows 2048 3 particles (~ 8 bil¬ 
lion) in a cubic box of side length 250 /i _1 Mpc using 
the ART code (Kravtsov et al. 1997). Its mass resolution 
(1.36 x 10 8 h~ 1 M Q ) and force resolution (1 h~ x kpc) 
allow it to resolve halos down to ~ 1O 1O M 0 . Bolshoi 
adopts a flat ACDM cosmology (Dm = 0.27, Da = 0.73, 
h = 0.7, ag — 0.82, n s m 0.95) which is very similar to the 
WMAP9 best-fit cosmology (Hinshaw et al. 2012). Halos 
were found using the Rockstar halo finder (Behroozi 


et al. 2013d), which associates particles to peaks in the 
phase-space density; the finder also makes use of tempo¬ 
ral information to improve stability in major mergers. 
Rockstar has been shown to have excellent recovery 
of both halo and subhalo properties (Knebe et al. 2011; 
Onions et al. 2012). Merger trees were constructed us¬ 
ing the Consistent Trees algorithm (Behroozi et al. 
2013c), which predicts the gravitational evolution of ha¬ 
los from one snapshot to another so as to identify and 
repair inconsistencies (see also Srisawat et al. 2013). 

Throughout the paper, we neglect subhalos and ex¬ 
clusively use host halos in all of our calculations. We will 
quote results for a halo’s mass accretion rate, dM v i r /dt. 
For a halo identified at cosmic time 1 1 , we compute this 
quantity as 

dMvir/dt = M ^( t l- T d y n(tl))-Mv ir (t 1 ) ; (1) 
Tdyn(tl) 

in Eq. 1, M v i r (ti — Td yn (ti)) is the virial mass of the 
halo’s main progenitor at time t\ — Td yn (ti) as identified 
by Rockstar, where Td yn = (GA V irPcrit) -1 ^ 2 , with A v i r 
being the virial overdensity (Bryan k Norman 1998). 


3 RESULTS 

In §3.1, we review the physics linking halo accretion rates 
to large-scale environment. We then present the primary 
result in §3.2: the mass accretion rates of pairs of ha¬ 
los are correlated for pair separations out to 10 Mpc, a 
phenomenon we dub halo accretion conformity. In §3.2 
we also demonstrate that the underlying physics of this 
effect is quite intuitive, and is due to the halo pair mutu¬ 
ally evolving in the same large-scale tidal environment. 
We show the redshift-dependence of halo accretion con¬ 
formity in §3.3. 

3.1 Tidal Forces and Mass Accretion Rates 

For a detailed review of the connection between the mass 
accretion rate of a halo and its large-scale environment, 
we refer the reader to Hahn et al. (2009). Here, we sum¬ 
marize the basic physical picture. 

Consider a pair of halos separated by distance D, 
calling the larger halo the primary and the smaller halo 
the secondary ; i.e., so that M pr i m > M sec . Because of the 
tidal force exerted by the primary upon the secondary, 
the only stable circular orbits around the secondary halo 
are for particles whose distance is smaller than the Hill 
Radius, 

RHill = D (Mgec/SMprim ) 1 / 3 

— Rsec(D/ V^LRprim) ( 2 ) 

The smaller the distance between the halos, the stronger 
is the tidal force, and the smaller is the radius around the 
secondary that circular orbits are stable (i.e., the smaller 
is Rmn). The Hill Radius is thus an upper bound on 
the spatial extent of newly infalling material that can 
remain gravitationally bound to the secondary halo. In 
this way, the presence of the primary regulates the mass 
accretion rate of the secondary. Figure 1 shows a cartoon 
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illustration of the Hill Radius for a pair of primary and 
secondary halos. 

For a secondary halo of virial radius R v ir, when any 
primary halo exerts a tidal force resulting in Rmw ~ R v ir, 
we should expect the secondary halo to experience a 
strong tidal field that inhibits its mass accretion rate. 
Motivated by this natural expectation, we quantify the 
strength of a halo’s large-scale tidal field with Rmn-min, 
defined as the minimum Rmw due to any of the halo’s 
more massive neighbors M > M sec , 

Rmil-min = MIN {i?Hill} M > M sec (^) 

In Figure 2, we show that the above natural expec¬ 
tation is borne out quantitatively. For z = 0 Bolshoi 
host halos in a narrow range of virial mass 2 M v i r ~ 
10 u Mq, we plot the mean present-day mass accre¬ 
tion rate dM v i r /dt of secondary halos as a function of 
Rmii-min/Rvir- The gray band in Figure 2 shows the 
error-on-the-mean, which we estimate here and through¬ 
out the paper via jackknifing 125 sub-volumes of the sim¬ 
ulated box. The blue shaded vertical region shows the in¬ 
ner 50% of the Rmii—min— distribution of 10 n M© halos. 

Figure 2 illustrates the following rough rule of 
thumb. For secondary halos with Rmn-min > 3R v ir, there 
is only a very weak dependence of dM v i r /dt on environ¬ 
ment. On the other hand, when Rmii-min < 3R v ir we can 
see that dM v i r /dt drops precipitously. On average, z — 0 
secondary halos with Rmn-min < R v ir have accreted no 
mass over the past dynamical time Td yn ~ 2Gyr. 

As shown in Hahn et al. (2009), there is a pro¬ 
nounced halo mass dependence to the trend shown in 
Fig. 2. For halos less massive than the collapse mass 
(Mcoii, which is ~ 10 12 ' 6 M© at 2 = 0), the mass accre¬ 
tion rate exhibits the sharp Rmn-min —dependence shown 
in Fig. 2; for more massive halos (M v i r > M co ii), there 
is only a very weak dependence. Relative to higher-mass 
halos, the accretion rates dM v i r /dt of lower-mass halos 
are more strongly affected by their large-scale environ¬ 
ment. We will see in the next section how this previously- 
established trend manifests in halo accretion conformity. 


3.2 Halo Accretion Conformity 

3.2.1 Definition & observational motivation 

As discussed in §1, recent SDSS measurements (Kauff- 
mann et al. 2013) have shown that the SSFR of a sample 
of “primary” central galaxies is correlated with the mean 
SSFR of all the “secondary” galaxies in the neighborhood 
of the primaries; the correlation persists out to projected 
separations of R p ~ 4 Mpc. Furthermore, the strength 
of this correlation weakens as the stellar mass of the pri¬ 
maries increases. 

In this section, we show that dark matter halos ex¬ 
hibit directly analogous trends in their mass accretion 
rates. Primary dark matter halos with above-average 
dMvir/dt tend to be located in an environment populated 
with secondary halos that are also rapidly accreting mass. 


2 All mass bins in this paper are 0.35 dex in total width. 
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Secondary Halo 

Msec 


Figure 1 . Diagram of the Hill Radius R mi i- Due to the 

tidal field produced by the primary, dark matter particles be¬ 
yond Rmw = D (M sec /SMprim ) 1 ^ cannot accrete onto and 
remain bound to the secondary halo. Thus halos in stronger 
tidal fields have smaller Hill Radii. When Rmii is roughly 
equal to R v i r , the virial radius of the secondary, we should ex¬ 
pect the secondary halo to experience a significant suppression 
of its mass accretion rate dM v i r /dt. 



Figure 2. Strong tidal fields stifle mass accretion rates. 

For each secondary halo with M v i r — 10 11 M©, we compute 
Rmil—min 5 the minimum Hill Radius due to any more massive 
halo in the simulation. The black curve shows the present- 
day mass accretion rate dM v i r /dt of the secondary halos as 
a function of Rmil-minj the gray band shows our jackknife 
estimation of the error-on-the-mean. The shaded blue region 
shows the inner 50% of the Rmil-min - distribution of 10 11 M© 
halos. The horizontal dashed line shows the median dM v i r /dt 
of all 10 11 Mq halos; the vertical dashed line guides the eye 
to the corresponding horizontal axis value. Tidal forces with 
Rhiil-min < 3R v ir dramatically stifle accretion onto the sec¬ 
ondary; for Rmil—min > 3R v ir i the large-scale environment 
does not significantly impact dM v i r /dt. 


The converse is also true for primary halos with below- 
average dM v i r /dt. We refer to this phenomenon as halo 
accretion conformity. 
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Under-dense Environment 


RHill-min from 

tertiary halo 






RHill-min from 
tertiary halo 



Over-dense Environment 



Figure 3. The physics of halo accretion conformity. The primary (red) and secondary (gray) halos are so distant that the 
direct tidal force between them has a negligible influence on either of their mass accretion rates. However, their accretion rates 
can still be correlated because the halos mutually evolve in the same large-scale environment. A halo evolving in a crowded field 
of other halos will tend to be located near another massive tertiary halo; this results in a strong tidal field, a small Ruilh and a 
suppressed mass accretion rate. Conversely, halos evolving in a void-like environment tend to be more isolated, with fewer massive 
tertiary halos available to exert a significant tidal force. Since this phenomenon pertains to both primaries and secondaries alike, 
there is a natural tendency for correlated assembly histories to arise between halos evolving in a similar tidal environment; this 
tendency is the physical origin of halo accretion conformity. 
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3.2.2 Basic physical picture 

Before presenting our measurements of conformity in sim¬ 
ulated halos, let’s first consider how the physics discussed 
in §3.1 should naturally lead to halo accretion confor¬ 
mity. Consider the cartoon diagram shown in Figure 3. 
We show a primary and secondary halo pair embedded in 
a large-scale environment of tertiary halos. The top panel 
illustrates an under-dense, void-like environment with a 
scarcity of tertiaries; the bottom panel an over-dense en¬ 
vironment crowded with many tertiaries. 

In the over-dense environment, both primary and 
secondary halos tend, on average, to be close enough to 
some other tertiary halo such that Rmii-min impinges 
on R v ir- Conversely, in the under-dense environment, the 
typical .Rh ill—min is substantially larger than R v i r for 
both primaries and secondaries. Coupling this observa¬ 
tion with Fig. 2, we can intuitively see that evolving in 
a crowded environment tends to stifle present-day accre¬ 
tion rates, and conversely for under-dense environments. 
This effect will result in primary halos being surrounded 
by secondaries with correlated dM v i r /dt, the definition of 
halo accretion conformity. 

Note that in this diagram, the tertiary halos produc¬ 
ing Rmii-min are distinct. As we will see, the dM v i r /dt 
correlation between neighboring halos is not due a direct 
tidal force on the primary and secondary by the same 
massive tertiary, nor due to direct tidal interaction be¬ 
tween the primary and secondary. We justify this inter¬ 
pretation of the physical cause of halo accretion confor¬ 
mity in §3.2.4. Before doing so, in the next section, §3.2.3, 
we quantify the scale- and mass-dependence of confor¬ 
mity in low-redshift halos. 

3.2.3 Conformity at low-redshift 

In Figure 4 we show the primary result of this work: 
the measurement of halo accretion conformity. Roughly 
speaking, Figure 4 shows the mass- and scale-dependence 
of how the mass accretion rates of lower-mass secondary 
halos “know about” the mass accretion rates of higher- 
mass primary halos. We separately study this signal for 
several different choices of primary and secondary halo 
mass, as shown in the different panels. In detail, we cal¬ 
culated the results in Figure 4 as follows. 

For each halo P[ in a given sample of primaries, 
we identify all secondary halos located in a spherical 
shell of radius D about the primary. We then calculate 
(dM v i r /dt)Pj (D), the mean dM v i r /dt of the secondaries 
in the shell that surrounds primary R. The vertical axis 
in all panels of Figure 4 shows the mean of all these ac¬ 
cretion rates, weighting each primary equally. We choose 
this method (rather than giving equal weight to each sec¬ 
ondary) so that primaries in over- and under-dense envi¬ 
ronments are treated on equal footing. We note, though, 
that our results are qualitatively insensitive to this choice. 

The black curves in Figure 4 show the results of the 
above calculation for all primary and secondary halos 
whose masses are indicated at the top of each panel. 
In direct analogy to the observational measurement of 
galactic conformity, we have repeated the above ex¬ 
ercise after first sub-dividing the sample of primaries 


into quartiles of primary dM v i r /dt. In each panel, the 
blue curve shows (dM v i r /dt) of secondaries surround¬ 
ing the fastest-accreting quartile of primaries; the red 
curve shows the same, but for secondaries surrounding 
the slowest-accreting quartile of primaries. 

The most important feature of Figure 4 is the sep¬ 
aration of the red and blue curves. This indicates that 
fast-accreting primary halos tend to be surrounded by 
an environment of fast-accreting secondaries, and con¬ 
versely for slow-accreting primaries. For lower mass ha¬ 
los, this separation persists to 10 Mpc, showing that halo 
mass accretion rates are correlated for spatial separations 
greater than thirty times the virial radius. 

In Figure 5 we illustrate the mass-dependence of halo 
accretion conformity in more detail. The vertical axis in 
Figure 5 shows the ratio of each blue/red curve in Fig¬ 
ure 4—i.e., a measure of how secondary halo accretion 
rates increase nearby fast-accreting primaries. Each curve 
in Figure 5 corresponds to a particular panel in Figure 4, 
as indicated by the legend. Larger values of the vertical 
axis in Figure 5 correspond to stronger halo accretion 
rate conformity. 

The mass dependence of conformity appears to be 
quite simple: lower-mass halos exhibit stronger confor¬ 
mity. This is to be expected based on the results in Hahn 
et al. (2009). Conformity is an effect caused by the large- 
scale environment, which has a proportionally greater ef¬ 
fect on halos of smaller mass. We note that both the mass- 
and scale-dependence of the signal shown in Figures 4 & 5 
closely resembles the trends seen in the SSFRs and gas 
fractions of SDSS galaxies (Kauffmann et al. 2013). 


3.2.4 Quantitative justification of physical picture 

While the physical picture sketched in §3.2.2 is intuitive, 
it is not obvious that the tidal forces typically arise from 
different tertiary halos. Indeed, primary and secondary 
halos that are both close to the same group or cluster 
may have suppressed dM v i r /dt due to direct tidal inter¬ 
action with the same massive halo. This proposed ex¬ 
planation differs qualitatively from the one illustrated in 
Figure 3; in this alternative scenario, massive groups and 
clusters play a privileged role in how they impact the as¬ 
sembly history of the low-mass halos in their ~ 5 Mpc 
environment. In this section we conclusively rule out this 
interpretation: the halos of massive groups and clusters 
do not have special significance in the physics driving halo 
accretion conformity. 

First, in Figure 6 we show the strong correlation be¬ 
tween large-scale environmental density and tidal field 
strength. For z — 0 Bolshoi halos in a 0.35 dex-width 
bin bracketing M v i r = lO n M 0 , we show the relation¬ 
ship between (Rmii-min) and environmental density. For 
our density estimator, we use 1 + ^ 5 , computed as the 
mass density in a 5 Mpc sphere divided by the cosmic 
mean matter density. As expected, the denser the envi¬ 
ronment, the smaller the typical value of Rh ill -min, and 
the stronger the tidal field. 

Second, in Figure 7 we show explicitly that a pri¬ 
mary halo’s tidal field is correlated with the tidal field 
of secondary halos out to large distances. We proceed in 
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Figure 4. Halo accretion conformity at z = 0. We compute dM®?f/dt of lower mass, secondary host halos in spherical shells 
surrounding higher mass, primary host halos. In each panel, the black curve shows the mean value of dM®^r/dt as a function of the 
distance between the primary and secondary. The masses of the primary/secondary pair are indicated by the title of each panel, 
where we have used mass-bins 0.35 dex in total width; gray bands denote the jackknife error-on-the-mean. Additionally, we have 
repeated each calculation after first sub-dividing the primary halos into quartiles of dM^f im /dt. The mean dM®? r c /dt of secondary 
halos surrounding the fastest- and slowest-accreting quartile of primaries is shown with the blue and red curves , respectively. The 
upshot is that i) fast-accreting primaries tend to be surrounded by fast-accreting secondary halos, and conversely; ii) both the 
mass- and scale-dependence of this signal closely resembles the conformity signal in SDSS observations of central galaxy SFR. 



Pair Distance [Mpc] 

Figure 5. M vir —dependence of halo accretion conformity strength. We plot the ratio of the blue-to-red curves in each of 
the panels in Fig. 4 , referring to this ratio as the “conformity strength.” Halo accretion rates are correlated out to R ~ lOMpc, a 
trend that weakens with increasing mass. 
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Figure 6. Denser environments have stronger tidal 
fields. For M v i r = 10 11 Mq halos at z = 0 , as a function 
of 5 -Mpc over-density we show the mean tidal field strength, 
as quantified by min- 



Figure 7. Tidal fields are correlated on large scales. 

We show the mean value of -Rhill—min of secondary halos in 
spherical shells surrounding primaries. In analogy to Fig. 4 , 
we have repeated the calculation after first dividing the pri¬ 
mary sample into quartiles of min - The mean Rhill—min 

of secondaries surrounding primaries in the top (bottom) quar¬ 
tiles of -^hill^min shown with the blue (red) curve. With 
the dashed, vertical black line we show the distance at which 
Rhill — 3R v i r (see Fig. 2). Primary halos in a strong tidal field 
are surrounded by secondaries that are also in a strong tidal 
field, the origin of halo accretion conformity. 

an exactly analogous fashion as we did to calculate the 
results shown in Figure 4. Here we focus on primary ha¬ 
los of mass M v ir = 10 12 Mq, and secondary halos of mass 
M v ir — 10 11 Mq. In spherical shells surrounding each pri¬ 
mary, we compute the mean value of R um -min of the sec¬ 
ondaries in each shell, and then average the shells. Thus 
the black curve in Figure 7 shows the average Rmn-min 
value of the secondaries as a function of their distance 
to the primary. The closer a lO n M0 secondary halo is 


to a 1O 12 M0 primary, the denser the environment of the 
secondary, and the stronger the tidal field exerted by the 
environment on the secondary. 

Again we divide our sample of primaries into quar¬ 
tiles, this time based on the Rmn-min value of the pri¬ 
mary. The mean tidal field strength of secondaries sur¬ 
rounding primaries in the weakest tidal field are shown 
with the blue curve ; results for secondaries surrounding 
primaries in the strongest quartile of tidal field strength 
are shown with the red curve. The separation of the blue 
and red curves in Figure 7 tells us that primary halos in 
a strong tidal field tend to be surrounded by secondary 
halos also experiencing strong tidal forces, an intuitive 
result. 

With the vertical dashed line in Figure 7, we show 
the distance at which Rmn = 3R sec • As shown in Fig¬ 
ure 2, this is the tidal field strength necessary to have an 
appreciable impact on dM v i r /dt. Thus in order for the 
halo pair to exhibit a direct tidal influence on each other 
that is sufficiently strong to impact dM v i r /dt, they must 
be separated by a smaller distance than what is shown by 
the vertical dashed line. The separation between the red 
and blue curves remains large for separations that vastly 
exceed this distance. It is then not possible that halo ac¬ 
cretion conformity is caused by direct tidal interaction 
between the halo pair. 

Finally, in Figure 8 we assess the plausibility of 
the alternative scenario described in the beginning of 
this section. Again we focus on secondary halos of mass 
M S ec = 1O 11 M0, where the effects are largest. First, with 
the black curve in the top panel , we show the mass func¬ 
tion of the tertiary halos that are responsible for the 
-Rmn-min of the secondaries; we have scaled this mass 
function by the abundance of the secondaries, so that the 
vertical axis in Figure 8 is in dimensionless units. The red 
curve shows the same result, but only for secondaries with 
-Rmii-min < 3R v i r ; we remind the reader that Figure 2 im¬ 
plies that this is the tidal field strength required to impact 
dMvir/dt. The similarity of the red and black curves in 
Figure 8 implies that the dM v i r /dt —stifling tidal forces 
experienced by the secondaries are not preferentially ex¬ 
erted by massive groups or clusters. 

In the bottom panel , we compute the fraction of pri¬ 
mary/secondary halo pairs for which the same tertiary 
halo is responsible for the Rmii-min of both the primary 
and secondary. Figure 8 shows how this fraction varies 
as a function of the distance between the pair. At dis¬ 
tances D « 5 Mpc, where both the tidal field strength 
(Figure 7) and mass accretion rates (Figure 4) are still 
significantly correlated, over 90% of primary-secondary 
halo pairs have distinct tertiary halos responsible for the 
most significant tidal field. This rules out the possibility 
that mutual proximity to the same massive object can ex¬ 
plain correlations between low-mass halo accretion rates. 
The singular presence of a massive group or cluster has 
no special significance in the physics of halo accretion 
conformity. 
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Figure 8. Massive groups and clusters do not have 
special significance in the physics of halo accretion 
conformity. Top Panel: For M sec = 10 11 M© secondary ha¬ 
los, the black curve shows the mass function of the tertiary 
halos responsible for Rum— min? normalized by the mass func¬ 
tion of the secondaries. The strikingly similar red curve shows 
the same quantity, but only for secondaries with Rum— min < 
3 -Rvir- Bottom Panel: We calculate the fraction of pairs with 
the same tertiary halo responsible for Rmil— min as a function 
of the distance between the halo pair. Taken together, the two 
panels imply that halo accretion conformity is not preferen¬ 
tially influenced by massive groups/clusters, nor is the signal 
due to mutual proximity to the same object, justifying the 
picture illustrated in Fig. 3 . 


3.3 Conformity at High Redshift 

We conclude our results by showing how halo accretion 
conformity evolves with redshift. We repeat the calcu¬ 
lation illustrated in Figures 4 & 5, but now including 
Bolshoi halos at redshifts z — 0.5,1 and 2. We focus ex¬ 
clusively on a fixed pair of primary and secondary halo 
masses of mass 10 12 M© and 10 n M©, respectively, where 
the z = 0 signal is strongest. These calculations allow us 
to estimate how the observed galactic conformity signal 
should evolve with redshift under the assumption that 
halo accretion rates are correlated with galaxy star for¬ 
mation rates (see §4.5 for further discussion). 

Figure 9 shows the direct analog of Figure 4, only 


showing high-redshift results. The strength of the signal 
substantially weakens as we go farther back in cosmic 
time. This weakening is more easily seen in Figure 10, 
which shows the high-redshift analog of Figure 5. 

This weakening is as expected. Since collapse mass 
Mcoii decreases with increasing redshift, by comparing 
the conformity signal for halos of fixed mass, at higher 
redshift we are studying halos with larger values of 
Mvir/Mcoii. We have already seen in Figure 5 that con¬ 
formity weakens for higher mass halos, a fact also seen in 
SDSS observations (Kauffmann et al. 2013). Moreover, it 
has already been shown in Hahn et al. (2009) that M co ii 
sets the natural mass scale for how large-scale environ¬ 
ment impacts halo assembly. It is therefore completely 
natural that the conformity signal weakens at higher red¬ 
shift for halos of the same mass. As discussed in §4.5, we 
fully expect the same to hold true for future measure¬ 
ments of galactic conformity at higher redshift, which 
would provide strong observational support for the natu¬ 
ral hypothesis that central galaxy SFR and dark matter 
halo accretion rate are correlated. 


4 DISCUSSION 

In §4.1 we describe how conformity provides informa¬ 
tion about the coupling between galaxy and halo growth; 
we also sketch how future measurements of galactic con¬ 
formity at both low- and high-redshift can be used to 
guide the modeling of feedback processes in hydrodynam- 
ical simulations and semi-analytic models. We outline the 
fundamental connection to assembly bias in §4.2. We de¬ 
scribe how our work suggests a connection between large- 
and small-scale conformity in §4.3, and discuss alterna¬ 
tive quantifications of our results in §4.4. We conclude in 
§4.5 by discussing predictions for galactic conformity at 
higher redshifts. 


4.1 The Coupling of Halo and Galaxy Growth 

In the same work presenting the discovery of large-scale 
galactic conformity in SDSS data (Kauffmann et al. 
2013), it was also shown that a current semi-analytical 
model (SAM) of galaxy formation (Guo et al. 2011) pre¬ 
dicted only a very weak conformity signal. Since our re¬ 
sults show that conformity exists in the accretion rates of 
dark matter halos, the failure of this prediction is curious. 
After all, the star-formation history of a central galaxy in 
a SAM traces the mass assembly history of its host halo. 

To understand this puzzling failed prediction, it is 
useful to look to a model which does predict strong levels 
of galactic conformity: age matching (Hearin & Watson 
2013; Hearin et al. 2014a; Watson et al. 2015). In the 
age matching model, older (earlier-forming) halos host 
central galaxies with older stellar populations relative 
to younger halos of the same mass. This naturally re¬ 
sults in galactic conformity: the present work shows that 
fast-accreting (slow-accreting) halos tend to congregate 
together, and halo ages and accretion rates are tightly 
anti-correlated at fixed mass (Hahn et al. 2009). 

Since the Hearin & Watson (2013) model presumes 
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z =1 z =2 




Figure 9. Halo accretion conformity at high-redshift. Same as Fig. 4, but for halos at different redshifts, labeled by the 
title of each panel. All curves pertain to primaries of mass M v i r = 1 O 12 M 0 and secondaries with M v i r = 1 O 11 M 0 at the indicated 
redshift. Halo accretion conformity generically weakens at higher redshift, due primarily to the time evolution of collapse mass. 



Pair Distance [Mpc] 


Figure 10. Redshift-dependence of halo accretion conformity strength. Same as Fig. 5, but for halos at different redshifts. 
Each curve shown in this figure is computed as the ratio of the corresponding blue-to-red curves in Figure 9. For redshifts z > 1, 
conformity on large scales (R > 3 Mpc) vanishes to undetectable levels. 
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that halo and galaxy age are in monotonic correspon¬ 
dence at fixed stellar mass, with no scatter, this model 
represents the limit of maximal coupling between the 
growth of a galaxy and its parent halo. In the opposite 
extreme limit, empirical models such as the Halo Occu¬ 
pation Distribution (HOD) and Conditional Luminosity 
Function (CLF) represent the limit of zero coupling, since 
in these models present-day virial mass alone governs the 
properties of the halo’s resident galaxies. HOD models, 
or indeed any model with zero coupling between central 
galaxy and halo growth, predict essentially zero two-halo 
conformity (Hearin et al. 2014b; Paranjape et al. 2015). 

Thus, the coupling strength between halo assembly 
history and galaxy star formation directly influences the 
predicted conformity effect. SAMs such as Guo et al. 
(2011) heavily rely on present-day halo mass (or gravita¬ 
tional potential) to determine star formation rates, as do 
many other current semi-analytical models (see Lu et al. 
2014a, for a review). The presence of history-unaware 
feedback will of course tend to erase the memory the 
galaxy has of its parent halo’s assembly, thereby eras¬ 
ing the conformity signal present in the underlying dark 
matter halo. 

We note that re accretion of ejected gas can also play 
an important role for galaxy formation at z — 0 (Nelson 
et al. 2015). If reaccretion were dominant, the same phys¬ 
ical mechanism—i.e., external tidal forces—would still 
cause galactic conformity because the tidal forces would 
affect gas reinfall times. In this case, the strength of con¬ 
formity would also depend on the turnaround radius of 
ejected gas. If the gas thermalized immediately or cy¬ 
cled nearby the galaxy, then external tidal forces would 
not play an important role, and no conformity would be 
present for the gas reinfall rate. However, if galactic foun¬ 
tains propelled gas to significant fractions of the virial 
radius, then stronger external tidal forces would result in 
gas taking significantly longer to reinfall, reducing the to¬ 
tal reinfall rate. Hence, in this scenario, conformity would 
be connected to the strength of star formation feedback. 

In either scenario, conformity represents a practically 
ideal statistic at z — 0 to constrain the coupling strength 
between galaxy environment and galaxy growth, as well 
as the channels modulating the strength of feedback. Ob¬ 
servational conformity measurements across cosmic time 
will provide invaluable insight into the physical processes 
regulating star formation and quiescence. 


4.2 The Fundamental Link between Halo 

Accretion Conformity and Assembly Bias 

Halo mass is the dominant variable determining the spa¬ 
tial distribution of dark matter halos. However, as dis¬ 
cussed in §1, the clustering of halos has an additional de¬ 
pendence on formation time Zf or m, a phenomenon dubbed 
assembly bias. This dependence can be quantified in 
terms of the relative strength of the two-point function 
of old and young halos of the same mass, as in Gao et al. 
(2005). 

This classical quantification of assembly bias is sim¬ 
ply linked to conformity. At fixed M v i r , earlier forming 
halos have lower dM v i r /dt at z = 0. Thus since early- 


forming halos cluster more strongly than late-forming 
halos of the same M v i r , then so do slow-accreting ha¬ 
los cluster more strongly than fast-accreting halos of the 
same M vir . Strongly clustered halos reside in preferen¬ 
tially denser environments, where tidal fields are stronger 
(Figure 6); the converse is true for weakly clustered halos. 
The basic result of this paper is that halos residing in the 
same tidal environment have correlated dM v i r /dt, and so 
we conclude that not only are conformity and assembly 
bias connected, they are alternative statistical quantifica¬ 
tions of the exact same phenomenon. 

4.3 Relevance to 1-Halo Conformity 

A correlation between the star-formation histories of 
neighboring galaxies was first discovered in Weinmann 
et al. (2006), who found that the SFR of satellite galaxies 
are correlated with the SFR of the associated central. By 
definition, central and satellite galaxies occupy the same 
dark matter halo, whereas the Kauffmann et al. (2013) 
measurements pertain to galaxies occupying distinct ha¬ 
los. Accordingly, Hearin et al. (2014b) dubbed these two 
signals “1-halo” and “2-halo” conformity, respectively. 

This paper is chiefly concerned with the large-scale 
2-halo conformity signal measured in Kauffmann et al. 
(2013). However, in light of the findings presented here, 
it is quite plausible that 1-halo and 2-halo conformity 
measurements probe the same underlying physics. Af¬ 
ter all, satellites lead most of their lives as centrals in 
the neighborhood of their ultimate host halo. Thus cen¬ 
trals living in host halos with below-average dM vir /dt 
will tend to accrete satellites that, prior to infall, have 
themselves already experienced below-average dM v i r /dt, 
and conversely. In fact, it is precisely this phenomenon 
which leads to the age matching model simultaneously 
exhibiting both 1-halo and 2-halo conformity (Campbell 
et al., in prep). In this way, 2-halo can naturally evolve 
into 1-halo conformity. 

4.4 Alternative Quantifications 

Throughout this work, we have focused exclusively on 
Rmii-min as the single number that encapsulates the in¬ 
fluence of a halo’s large-scale environment on its mass 
assembly. In developing our results, we explored a broad 
range of alternatives that could act as regulators of mass 
accretion. Our investigation of these alternatives included 

(i) distance to the nearest halo of (variable) mass M v i r ; 

(ii) large-scale density (with a variable kernel); 

(iii) host halo concentration, Vmax, and Zf orm ’, 

(iv) whether or not the host halo has been previously 
ejected from another halo. 

Of these alternatives, we find that Rmii-min is the 
strongest predictor of halo accretion rate, regardless of 
the summary statistic of the dM v i r /dt distribution (me¬ 
dian, mode, size of the low-dM v i r /dt tail, etc.). In addi¬ 
tion to the clear physical motivation for Rmii-min, this 
proxy has the distinct advantage of being computable in 
a galaxy catalog, at least in principle. Redshift-space dis¬ 
tortions will introduce noise into using Rmn-min as an 
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environmental proxy, but otherwise the only required in¬ 
gredient is a model for the M* — M v i r relation. Existing 
models (e.g, Behroozi et al. 2013a; Kravtsov 2013) may 
be sufficient for this purpose, since the mass-dependence 
of -Rmii-min is weak (see Eq. 2). 

We have also explored the dependence of halo accre¬ 
tion conformity on the timescale r over which dM v i r /dt 
is defined (see Eq. 1). We find that the signal decreases 
in strength for r < 1 Gyr, and essentially vanishes if the 
“instantaneous”, snapshot-to-snapshot timescale is used. 
Since conformity in galaxy SFR is realized in the real Uni¬ 
verse, this observation could signal a natural timescale 
over which galaxy and halo growth are correlated, or it 
could simply be a time- and/or mass-resolution effect. We 
leave this investigation as a subject for future work. 

4.5 Predictions for Conformity at Higher 
Redshifts 

The qualitative redshift-dependence of conformity in 
Figs. 9 and 10 translates straightforwardly to galaxies: 2- 
halo galactic conformity should be much weaker at higher 
redshifts. The logic of this expectation is quite simple. 
First, the stellar-to-halo mass relation evolves only very 
weakly with redshift (e.g., Behroozi et al. 2013b). Thus 
by studying the redshift evolution of conformity at fixed 
M v ir, we are effectively studying the redshift-evolution at 
fixed stellar mass. Second, the intrinsic mass scale in the 
problem is M co ii, the halo model collapse mass, which 
decreases monotonically with increasing redshift, so that 

M v ir M v ir 

M co n(z > 0) > M co n(z 0)' 

Finally, the influence of environment on a halo generi- 
cally decreases as M v ir/M C oii increases. Therefore, if the 
SSFR of a central galaxy is correlated with its host halo 
accretion rate, the prediction that conformity weakens at 
higher redshift follows inexorably. 

The largest halo accretion conformity signal in this 
study occurs for 10 12 M© primary halos with 10 11 M© sec¬ 
ondary halos . 3 These correspond to ~ 10 10 — 10 lo ' 5 M© 
primary galaxies with 10 9 ‘ 5 — 1 O 1 O M 0 secondary neigh¬ 
bor galaxies (roughly independently of redshift for z < 4; 
Behroozi et al. 2013b). Galaxies of these masses are cur¬ 
rently probed with precise redshifts out to z ~ 0.5 in, e.g., 
the PRIMUS (Moustakas et al. 2013), VIPERS (Garilli 
et al. 2014), and DEEP2 (Bundy et al. 2006) surveys. 
Hence our high-redshift predictions can be tested even 
with present datasets. 

If the strength of galactic conformity is correlated 
with the strength of halo accretion conformity, then Fig. 
10 provides a basic guide for how to relate z = 0 con¬ 
formity signals to z > 0 conformity signals. For example, 
Fig. 10 indicates that the conformity strength observed at 
a separation of 4Mpc at z=0 should be matched by a sig¬ 
nal of the same strength at a separation of ^2.2 Mpc (co¬ 
moving) at z = 0.5. The following fit well-approximates 

3 We note that we expect halo accretion conformity to con¬ 
tinue to increase in strength towards lower masses, but this 
remains to be seen with a higher resolution simulation. 


the expectation from Fig. 10: 


S g (R,z) 

« l + A(z)(R+1.65) B(z) 

(4) 

A(z) 

= 15.8 - 6.5z 

(5) 

B(z ) 

= -2.8 - 0.67* 

( 6 ) 


where R is the comoving separation between primary and 
secondary, z is the redshift of the observation, and S g is 
the strength of galactic conformity, measured in the same 
fashion in Figures 5 V 10, i.e., the upper-to-lower quartile 
ratio. 

We note that the prediction given in Eq. 4 is only 
approximate, particularly because we have only stud¬ 
ied conformity in three dimensional shells, rather than 
2 + 1-dimensional cylinders in redshift-space. Redshift- 
space distortions will generically reduce the signal, so the 
overall magnitude of the signal shown in Figures 5 V 10 
represents an upper bound. However, we have not incor¬ 
porated this effect into Eq. 4, as the appropriate redshift 
window depends on the survey used. Since the present 
literature does not even contain an order-of-magnitude 
estimate on how conformity should scale with redshift, 
we consider Eq. 4 to be a useful guideline until more de¬ 
tailed models and measurements become available. 

For primary galaxies with stellar mass larger than 
10 1 X M® (hosted by ~ 10 13 M© halos), Fig. 5 suggests 
that the conformity signal is roughly half of that for the 
lower mass ( 10 10 — 1 O 1 O ‘ 5 M 0 ) primary galaxies consid¬ 
ered above. This would suggest that the 2-halo galac¬ 
tic conformity signal would be marginally detectable at 
z m 0.5, and only at closer distances to the primary galax¬ 
ies (~1 Mpc). However, the virial radius for a 1 O 13 M 0 
host halo is already ~ 0.6 Mpc, and because galaxies 
continuously pass in and out of the virial radius, large ha¬ 
los can directly affect quenched fractions at several times 
their virial radii (Wetzel et al. 2014). Hence, it is likely 
that any signal found at higher redshift would be hard 
to attribute solely to 2 -halo galactic conformity (see also 
Paranjape et al. 2015, for elaboration of this point). Fi¬ 
nally, we note that if 1 -halo and 2 -halo conformity are 
indeed manifestations of the same underlying physics, as 
discussed in §4.3, a corollary to our findings is the predic¬ 
tion that 1 -halo conformity should also vanish for z > 1. 


5 CONCLUSIONS 

Our primary findings are: 

(i) The mass accretion rates of dark matter halos are 
highly correlated on large scales, as shown in Figures 4 & 
5. We dub this phenomenon halo accretion conformity. 
Both the scale- and mass-dependence are qualitatively 
similar to SDSS measurements of galactic conformity, 
strongly suggesting that at fixed mass, dark matter halo 
accretion rate and central galaxy SFR are correlated. 

(ii) The origin of halo accretion conformity is the mu¬ 
tual evolution of halos in the same large-scale tidal field 
(see Figures 6 V 7). Massive groups and clusters play no 
special role in the physics giving rise to the signal (Fig. 8 ). 

(iii) We make quantitative predictions for the time 
evolution of galactic conformity (Eq. 4), suggesting that 
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the signal will drop rapidly at higher redshifts. Our pre¬ 
dictions are testable with existing datasets to z ~ 1. Fu¬ 
ture conformity measurements will be highly informative 
about the physical processes driving galaxy evolution. 
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